On exact solutions of Maxwell-Bloch system 
for two-level medium with degeneracy 



S. B. Leble*, 
Wydzial Fizyki Technicznej i Matematyki Stosowanej 
Politechnika Gdan'ska, ul. G.Narutowicza, 11/12 80-952, Gdan'sk-Wrzeszcz, Poland, 

email leble@mifgate.pg.gda.pl 
N.V.Ustinov 

Kaliningrad State University, Theoretical Physics Department, 
Al.Nevsky st.,14, 236041 Kaliningrad 
email ustinov.leble@theor.phys.ksu.kern.ru 
*permanent adress 



Abstract 

Maxwell-Bloch system describing the resonant propagation of electromagnetic pulses 
in both two-level media with degeneracy in angle moment projection and three-level 
media with equal oscillator forces is considered. The inhomogeneous broadening of 
energy levels is accounted. Binary Darboux Transformation generating the solutions 
of the system is constructed. Pulses corresponding to the transition between levels 
with largest population difference are shown to be stable. The solution describing the 
propargation of pulses in the medium exited by the periodic wave is obtained. The 
hierarchy of infinitesimal symmetries is obtained by means of Darboux transformation. 



1 Introduction 



Investigations of laser emission interaction with resonant media sufficiently promoted the 
understanding and mathematical description of this phenomena 1 . The important kick for it 
gave the invention of the self-induced transparency (SIT) phenomena 2 that is a passage of 
a powerful and ultrashort pulse without loss of a form and energy at a resonant medium. 
The group velocity of such pulses depends on their amplitude and duration: shorter pulse 
has larger velocity. It was found that while the propagation of pulses of different duration 
there may be situation when a quick pulse reaches more slow one, after "collision" both 
go out without change of their form and velocities only phase shifts appear. From mathe- 
matical point of view this wonderful stability of SIT pulses is the corollary of the complete 
integrability of Maxwell-Bloch (MB) equations that describe SIT phenomenon by two-level 
medium approximation with nondegenerate levels 3 , recently developed for the case of two- 
photon processes 4 . The SIT pulses correspond to soliton solutions of the equations and their 
collision properties reflect the asymptotic properties of multi- soliton states with respect to 
one-soliton ones. 

For real media the energy levels are usually degenerate that contribute to peculiarities of 
polarizied pulses propagation. In the Ref. 5 the integrability of MB equations for arbitrary 
polarization of light being resonant with quantum transitions % — — > j a = 1 have been 
established. It was shown that solitons with circular polarization are stable. 

This problems of interest relates to the interaction of many-frequency laser pulses with 
a resonant many-level medium. For example such many-frequency action widen possibilities 
of isotop dividing or simulation of chemical reactions as well as spectroscopy investigations. 
The peculiarities of coherent pulses that propagate in many-level media may be used in the 
frequency transformation effectivisation 6 . 

Let's consider a medium that may be described in the two-level atoms approximation 
with energy levels E a ^ that are degenerate by projections m and ji of total momenta j a ,b- 
Let plane electromagnetic wave propagates along z and its electric field has the form 



where t is time. The carrier frequency uj = kc of the pulse is close to u Q = 2n(E b — E a )/h of 
the atom transition j a — > jt>. 

The evolution of the pulse (1) in the semiclassical approach and coherent approximation 
is described by generalized MB equations 5 that may be represented in the dimensionless form 

as 



E = 'E exp[i(kz — out)] + c.c. 
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In those equations the amplitude E Q of a light pulse, coefficients of the density matrix 
Pw, are connected with the functions e q , R^y, R m ,m>, R^m- 

E q = he q /2nt d, 
Pw' = N f(v)R^', 

Pm,m' -N f{jf) Rm,m' i 

P^ m = N f(r])R^ m exp[i(kz - cut)] 

here d is the reduced dipole moment of the transition j b — > j a , N Q - density of resonant 
particles, t = (3h/2nud 2 N ) 1 / 2 is time dimension constant. A = {u — cu )t - resonance 
shift, r = (t — z/c)/t ,( = z/ct , 



Jim = (-I)' 6 "" 
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by the index q the spherical components of the vector e are labelled . 
Let's consider the transition j b = 0— > j a = 1 and introduce 

e + = — ie + i,e_ = — ie_i, 

IT-b Rfxtxi \n=ix'=0i 
11(1+ Rmm'\m=m'=lj Tla— Rmm'\m=m' = — 1 
Rmm'\m=m'=—li ^+ Rfj/m I /i=0,m=l j ^— Rfj/m I [i=0,m=— 1 (^•^) 

The quantization axis is choosen along the light propagation direction (axis () 

As the interaction of a quantum system with the transverse electromagnetic waves is 
accompanied by transitions with angular momentum change (-1 or +1), the MB system may 
be simplified. 

e_ c = - < v- > 

e+,c = -<*/+> (1.3a) 
n a ~ )T = —v~e*_ — v*_ e_ 
n a +,r = ~v+e\ - v* + e+ 
n b,r = ^+ e + _ u X e + + f-e* + z/*e_ (1.36) 

v-,r = ~ + ( n a~ - n b )e- + v a e + 

v+ )T = -i(r) ~ A)z/+ + {n a+ - n b )e + + z/ a e_ 

v a ,r = -v+e*_ - v*_e+ (1, 3c) 

Here the angle brackets denote the inhomogeneous broadening averaging with normalized 
partition function f(rj) 

/oo 
f(v)Fdrj . 
-oo 



If the resonant atoms ensemble is in pure state, one may rest only five equations in the 
system (1.3). 

= — < a 3 a{ > 

= — < a 3 a* 2 > (1.4a) 

Oi,t — ~~ A)ai/2 — a 3 e*_ 

«2,r = i(v - A)a 2 /2 - a 3 e* + 

a 3,r = i(v - A)a 3 /2 + a x e- + a 2 e + (1-4&) 

where 01,2,3 are amplitudes of the probability density (population numbers) of the lower 
and two upper levels. 

The same system of nonlinear equations describes the resonant propagation of two- 
frequency radiation with equal shifts from resonance in three-level medium where one of 
the transitions is forbidden and the oscillator forces of the rest two transitions are equal. 



Zero curvature representation, reductions of poten- 
tials and Darboux transformation 



The Maxwell - Bloch equations (1.2) represents a complicated system that is not solved 
analytically in the case of arbitrary angular momenta. However in the case of resonant 
transitions % — > j a — 1; % — > j a — the equations (1.2) may be presented as the compatibility 
condition of some linear system (zero curvature representation) with a sperctral parameter 
(SP) A. For the case of equations (1.3) it is the following system 

^ =< a(X)A > ip (2.1) 
where V — U — XJ,J — diag{l, 1, -1}, a(X) = (2A + i(rj - A))" 1 , 



/ -e* 
U= -e* + 
\ e_ e + 




From the obvious equality ^ T = ip^ T we obtain the system 

U c = [J,<A>]/2, 

A T = [U,A]-i(r j -A)[A,J]/2, (2.2) 

from which it follows that the elements of U and A matrices should satisfy MB system (1.3). 
If matrix elements of A are representable in the form of Akj = ak<i*(k,j = 1, 2, 3) (the pure 
quantum states) then from (1.3) the system (1.4) follows. 



The equations (1.3,4) appear as the compatibility conditions of the conjugate Lax pair as 
well 

£ f = -f < a{n)A > (2.3) 

where W = U — kJ. 

As the matrices U and A satisfy the following reduction 

{ U A + S:° , m 

one may check that there in the space of matrix solutions of ZS equations exists the autho- 
morphism in the sense of the Ref. 7, 

#(A) -> [*(-A*) + ] _1 e *(A). (2.4) 

By the way there is the coupling between the WF spaces of right and conjugate (left) zero 
curvature representation 

K=-\*,i = ^ + (2.5) 

Such functions are refered as coupled by an authomorphism. 

The following statement is valid: The equations (1.3) are the corollary of (2.2) when the 
additional condition (2.4) is posed. 

The DT technique that is developed in the Refs. 8,9 may be applied for the construction 
of solutions hierarchy for the MB system (1.3). We specialize the scheme here. Let and 
X are functions of the right and left (conjugate) problems (2.1,3). Spectral parameters (SP) 
are \i and v correspondingly. We also suppose that those functions are connected by the 
authomorphism and 

v = -ff,X = </>*• 
The transformed WFs by the binary 8 DT 

^l] = [l-(/x + //)P/(A + //)]V (2.6a) 

£[1] = £[1 - (/i + mu*)P/(K - /!)] (2.66) 

are also the solutions of right and conjugate ZS zero curvature reprezentations with new 
potentials V[l] = U[l] - AJ, W[l\ = U[l\ - kJ and 

A[l] = A - 2(/i + //)[«(//) AP + a(fi)*PA - (a(fi)A + a(fi)*)PAP], (2.6c) 

where 

U[l]=U-( f jL + l jL*)[J,P]. (2.6d) 
The matrix P at the equations (2.6) is defined in analogy with Ref. 8 by 

Pkj = 0fe0*/(0 + ? 0) 

As P + = P the reduction (2.4) is conserved. It means that elements of matrices U[l] and 
A[l] give new solutions of MB equations (1.3). For the transformed potential coefficients that 
have the sense of dimensionless electromagnetic field components the following expressions 
are valid 

e_[l] = e_ + 2(/i + /i*)030!/(0 + , 0) (2.7a) 



e+[l] = e + + 2(» + /i*)0 3 0;/(0 + , 0) (2.76) 

If the initial state of the quantum system was pure the matrix Akj = cika* and the result of 
the transformation is pure as well with 

a[l] = a - 2(/i + /j>(/i)*Pa, 

the normalization is conserved too. 

As by the map (2.6) the transforms of WFs coupled by an authomorphism should also be 
connected by £[1] = ip[l] + the reduction (2.4) will be valid by iterations of BDT. Final result 
may be expressed via 0^ and that are WFs of initial zero curvature representation to 
be coupled by authomorphisms in pairs 9 . 

3 Solutions over the zero and periodic backgrounds. 

In this section the construction of solutions of MB (and NS) systems by the transformation 
(2.6,7) will be considered. 

3.1 Starting from zero functions e_ = e + = 0,n a+i _ — v a — z/ +i _ = 0, but rib = 1 
we obtain after first iteration the one-soliton solution that generalize 2ir pulse of two-state 
system. This solutions have been built in the paper 5 within the 1ST method. However a 
construction of two-soliton solution met difficulties and the analysis of soliton collisions was 
made by asympotic methods 6 . By the way in the Ref. 10 the algebraic approach that have 
features of Backlund and "Dressing" methods was developed and the breather two-soliton 
(with equal velocities) solution have been found. 

After a second iteration preserving the reduction (the WF have SP coupled by ix^) — A*(i)*) 
we obtain the solution of MB system with potential elements 

e_[2] = — 2(aiCiA 2 exp[— irj\ — a 2 c 2 Ai exp^r?] — aic 2 A* — a 2 c^A)/A[2] 
e+[2] = -2(b lC \A 2 exp[-^] - fccJAi expH - b lC * 2 A* - 6 2 c*A)/A[2] 

A[2] = AiA 2 - | A | 2 
Ai = [(| a, | 2 + | h | 2 )exp[-tf]+ | Cl | 2 exp[tf]]/(2/^) 
A 2 = [(| a 2 | 2 + | b 2 | 2 )exp[-^]+ | c 2 | 2 exp[tf]]/(2/^) 
A = [(aia 2 + bib* 2 ) exp[-^ - itheta] + c x c\ exp[^ + i9]]/(2ix) 

where 

■& = 2fji R (T+ <\ a{n) | 2 >£), 

6 = 2iijT- < (2fn + 7] - A) | a(n) | 2 > 0, 

where ix — /ir + i/ii,ai,2,bl,2,ci >2 are complex constants. The transformation of such a 
structure we may also name the binary one 1 ' 8 . This solution generalize the solution from 
Ref. 10. It is characterized by eight real parameters: ixr, /ij, distance between pulse centers, 
electric field components ratio of the pulses and phase shifts of the same component of both 
pulses. At Ref. 10 the solution contain four parameters with fixed amplitude and phase shift 
ratios. 



3.2 Let us consider the case of arbitrary level populations ("nonzero dipole temperature" 
of the medium): n a +n a ^n b ^ 0, e_ = e + = V- = u + = u a = 0. WFs of the sistem (2.1) have 
the form 

■01 = Ciexp(— Ar+ < a(A)n a „ > C) 
•02 = C 2 ea;p(-Ar+ < a(X)n a+ > () 
"03 = C 3 exp(\r+ < a(X)n b > () 
Transforming by (2.6) with WF and spectral parameter fi one get 

e_[l]= W30!/(0 + ,0), 

e+[l]= W 3 0;/(0 + ,0), 

where 

01 = Ci exp(-/ir+ < a(//)(n a _ - n a+ ) > C/2) 

02 = C 2 exp(-/ir+ < a(^)(n 0+ - n a _) > C/2) 

03 = C 3 exp(/rr+ < a(ii)(2n b - n a _ - n a+ ) > C/2) 

Let rife > n a+ > n a _. Then while r — > — oo,e+[l] — > at arbitrary C and if r — > +oo, 
then e_[l] — > 0. So the solution describes a transformation of a pulse resonant with the 
transition which has the minor population difference to one that corresponds to the larger 
population difference. It means that such pulses are not stable. In the paper 6 that process 
the more cumbersome 1ST technique in terms of Riemann - Hilbert problem we see the 
opposite statement . Our result seems to be more physical: the transition from the basic 
state to more populated one would be more transparent. 

3.3 Now we shall construct solutions of the MB system starting from a periodic back- 
ground. Let the plane electromagnetic wave is propagated along one of transitions 

e + = Eexpz(k( + wr), e_ = 0, 

where E is the a complex constant. As in MB equations the shift from resonance is taken 
into account we may put uu — 0. Let also the second level is empty n a _ = 0,nt + n a+ = 1. 
Then from the equations (1.3) one have 

v- = v a = 

v + = t(2n b — 1){uj + i] — A) _1 Eexp«(A;C + wr), 

k — — < (2n b -l)(u + rj- A) -1 >, 

n b = (l+(u + r ] -A)(4\E\ 2 +(u + r ] - A)' 1/2 )/2. 

During the derivation of the lower level population n b it was supposed that the quantum 
system had gone in the considered state from one with free upper levels. This condition 
fix the trace of the matrix A 2 that does not depend on r. It may be seen from the second 
equation of the system (2.2). The case of the arbitrary higher level population may be 
considered analoguesly. Relaxation terms from spontaneous transitions may be neglected if 
the dimensionless electric field amplitude | E \ is big enough. The solutions of the linear 
system are: 



■01 = C\exp{— At), 



■02 = (C + exp(-&) + C^exp(--&))exp((< a(X) > -ik)(/2) 
^ 3 = -(C+(A + <r))exp(7?) + C_(A - a)e:rp(-tf))e;rp((< a(A) > +ik)(/2)/E*, 

where 

= cr(r + i < (2n b - l)a(A)(?7 - A)" 1 > C) 
o- = (A 2 - | E \ 2 yi/2. 

Transforming by (2.6) and introducing a new parameter 7 = 7i? + iji instead of the SP \i by 

cosh(7) — fj,/ \ E \ 

one obtains 

e_[l] = -4£cosh( 7R ) cos(7/)030i exp(i/</2)/(0+ 0), 
e+[l] = -£(1 -4 J Bcos/i(7 R )cos(7 / )0 3 2 )ea:p(iA;C)/(0 + 0), 

where 

0i = Ci exp(— I £ I cosh(7 R ) cos(7/)(r+ < D > () + 
i(\ E I sinh( 7fi ) sin( 7j )r)- < D(r] - A) > C/2)). 
2 = C+ exp(^) + C. exp (-7?). 
03 = -£(C+ exp(^ + 7) + C_ exp(-^ - 7))/ I E |, 
and = i?ei? + ilmfi, 

Red =\ E \ cos(7/)(sinh(7 jR )r+ 

< (2n b — 1)(77 — A)" 1 ((t7 - A) sinh( 7R ) — 2\E\ sm( 7z )) J D > C), 

/mi? =| E I cosh(7 R )(sin(7/)r+ < {2n h - 1) (77 - A)~ 1 ((?] - A) sin(7/) + 

+2 I £ I sinh( 7i i))D > C), 

D = 2 I E | 2 (cosh(27 R ) + cos(2 7/ )) + (77 - A) 2 

This solution is parametrized by four complex parameters E, 7, C 1: C+/C- . The obtained 
solutions describe proceses of fission and fusion of radiation pulses that are resonant with 
different transitions. It shows the existence of pulses with different propagation velocities over 
the periodic wave background. The propagation of SIT pulses at arbitrary level poulation 
also demonstrate processes of frequency transformation. In order to verify this it is necessary 
to consider real parts of exponents indices in the expressions for 1? 2 ,03. The parameter 
7 is introduced namely for a convenience of this operation. The expressions that determine 
the character of pulses behaviours are however complicated but allow to plot electric field 
components. It is seen that the amplitude e + tends to | E \ at large ( and the phase of the 
periodic wave undergo a shift. 

The technique developed here may be applied for many-level systems in the cases when 
the corresponding MB equations are integrable. The studying of more complicated reductions 
of the ZS problem potencials seems to be interesting looking for more weak constraints on 
quantum systems parameters. Another possibility for applications is investigation of the 
conditions of pulse propagation in thin films as well as relection from them. 



4 Infinitesimal parameters of DT. 



In this section the solution for a linearization of the equations (2.2) to be obtained . The 
dependence of this solution on an arbitrary function give the possibility to widen classes 
of solutions of initial - boundary problems, looking for sequence of infinitesimal symmetries 
11 and for the stability of soliton solutions analisys. The main observation that allows to 
develop the technique, is that the nontrivial binary transformations (2.6) may generate a 
new potential that coincide with the seed one 12 . In this article BDT (not Matveev one) 
and the method of keeping reduction were used. It means that the DT is infinitesimal if 
U[l] = U, if v = fi and (x, <p) + 0. 

Linearizing the system (1.7) with respect to perturbations JjW and in relations to 
solutions U and A we get 



AW 



-[J,<AU>] 

^(77 - A) J + U , 1 + [£/ (1) ,A' 



(4.1) 



We shall suppose that v — /j, + 5, and 5 is small . 

Expanding (2.6) in Taylor series by the parameter 5 we get 

r u[i] = u + 8[j,p<®] + ... 

\ A[l]=A + 26a(fjL)[A,P^] + ... ' 

where P^ — ip <g) X j V 9 an d X are solutions of the direct and conjugate problems (2.8,9) with 
the spectral parameter /x, such that (x, if) = 1. 

From these results the following formulas for small disturbances are derived 



\ AW =2a(n)[A, P(°) 
After iterations of binary DT with the conditions 

jyil) = ^1) + §(1) , 



(4.2) 



where 
and 
one get 



<*(«)= O (tfM) (g = l,iV) 
« , 



jv 



JV 



,(0 V (0 



1=1 



where /3j — are constants. 

Going from sums to integrals we obtain the following representations of the solutions of 
the system (2.2) 



A w 



J. 



#(A)$(A)ciA 



Jc J 
A , f a(A)^(A)$(A)dA 



(4.3) 



where ^(A) and $(A) — are matrix solutions of the basic and conjugate problems (2.1,3) 
with a spectral parameter A. 

The validity of the formulas (4,2,3) one may check by the direct substitution in the 
equation (4.1). 

We notice also that the reduction (2.5) is compatible with the equation (3.3). By means 
of equality (3.5) it is possible to obtain solutions that satisfy the reduction, if integrals in the 
right-hand-side would be presented as a sum by contour integrals that may be parametrized 
in the following way 

/ *(A)$(A)dA = f 2 *(g(t))*(g(t))g'(t)dt+ 

+ [ t2 <S> + (g(t))* + (g(t))g*'(t)dt, 
Jti 

[ a(X)^(XMX)d\= f 2 a(g(t)Mg(t))<S>(g(t))gi(t)dt- 

- f 2 a*(g(t))<S> + (g(t))* + (g(t))g*'(t)dt . 
Jti 

The formulas of this paragraph are novel ones. As these expressions for small pertur- 
bations contain arbitrary matrix solutions of both problems (2.1,3) they may be applied 
for studying of stability of soliton solutions of the MB equation. It is also possible to con- 
struct infinitesimal symmetries 11,12 of the equations . and widen classes of initial - boundary 
problems solutions as well. 

5 Conclusion 

We would repeat that in 12 the infinitesimal binary DT was introduced and exploited to 
extract hierarchies of symmetries The technique of DT with reductions described above may 
also be applied to the case of NSS (Manakov equation) 13 as the first equation and the reduc- 
tion of the potential are the same. The only difference is in formulas (2.6a,b) and (2.7): the 
vector WF should be replaced by solutions of the consistent Lax problem 12 . Analoguesly one 
can built solutions of the NSS over the periodic background (the soliton case is considered 
at Ref. 6). Let 

e_ = 0, e + = E exp(i(k( + wr)) 

Then from Lax equations it follows that k = uo 2 — \ E | 2 . Solving the Lax system one get 
expressions for WF 

ipi = Ciexp(—Xr + iX 2 (), 
i> 2 = (C+exp(<d) + C-exp(-$))exp(-i(k( + cjt0)/2) 
^ 3 = -(C+(A + a - ik/2))exp($) + C_(A -a- ik/2)exp(-<d))exp(i(k( + cut), 

where 

$ = a (r + i(X + ik/2)C) 

a = ((A - ik/2) 2 - | E | 2 ri/2. 

If you put here A = ji and substitute the obtained formula for in the equation (2.7) you'll 
get the desired solution. Its difference from the solution of MB system is only in other 
dependence on the variable (. 



It is possible to obtain the integrable equation that combine systems MB and NS. This 
new integrable system may describe a propagation of polarized waves in Kerr nonlinear media 
with resonant atoms admixture similarily to linear polarized case of Ref. 14. The similar 
approach to perturbation theory based on Riemann - Hilbert problem is presented at the 
recent paper 15 . 
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